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Abstract. By using a generalization of Sturm-Liouville problems in discrete 
spaces, a basic class of symmetric orthogonal polynomials of a discrete variable 
with four free parameters, which generalizes all classical discrete symmetric 
orthogonal polynomials, is introduced. The standard properties of these poly- 
nomials, such as a second order difference equation, an explicit form for the 
polynomials, a three term recurrence relation and an orthogonality relation are 
presented. It is shown that two hypergeometric orthogonal sequences with 20 
different weight functions can be extracted from this class. Moreover, moments 
corresponding to these weight functions can be explicitly computed. Finally, 
a particular example containing all classical discrete symmetric orthogonal 
polynomials is studied in detail. 



1. Introduction 

Some special functions of mathematical physics such as classical orthogonal poly- 
nomials and cylindrical functions [29], are solutions of a differential equation of 
hypergeometric type [28l [29l [32] 

(1) 'J{x)y"ix) + T{x)y'{x) + Xyix) = 0, 

and extendible by changing equation ([!]) to a difference equation of the form 



(2) a{xis)) 



A 



Va;i(s) 



Vy(s) 



Va;(s) 



f(x(s)) 



Ay{s) , Vy{s) 



where 

Ax{s) = x{s + 1) - x{s), Va;(s) = Ax-(s 



Ax{s) Vx{s) 



A 



Xy{s) = 0, 
fis + l)-m 



Ax(s) '' ^ ' x{s + l) — x{s)' 

a{x{s)) and f{x{s)) are polynomials of degree at most two and one, respectively, 
in x{s), A is a constant, and xi{s) = x(s + 1/2). 
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The difference equation ([2]) , which is obtained by approximating the differential 
equation ([T]) on a non-uniform lattice, is of much importance |32| as its partic- 
ular solutions have been applied in quantum mechanics, theory of group repre- 
sentations and especially computational mathematics, where one can point to the 
Clebsch-Gordan and Racah coefficients with wide applications in atomic and nu- 
clear spectroscopy. There exist different approaches for the analysis of orthogonal 
polynomials of a discrete variable running from the classical references [HI I17j to 
the recent monograph |16| . which is a basic reference on orthogonal polynomials. 

Also there exists a number of numerical and symbolic methods for solving hy- 
pergeometric equations of type ([T]) or ([2]), which are of interest in applications, par- 
ticularly for cases containing symmetric solutions, such as resolution of the Gibbs 
phenomenon |10| lllj , Fourier-Kravchuk transform used in Optics [6] , approxima- 
tion of harmonic oscillator wave functions |4] , tissue segmentation of human brain 
MRI through preprocessing [T], reconstructions for electromagnetic waves in the 
presence of a metal nanoparticle |27| , efficient determination of the critical parame- 
ters and the statistical quantities for Klein-Gordon and sine-Gordon equations with 
a singular potential [7] , image representation [33l [M] and quantitative theory for 
the lateral momentum distribution after strong- field ionization [9j . 

The main aim of this paper is to introduce a basic class of symmetric orthogonal 
polynomials of a discrete variable with four free parameters, which is the polyno- 
mial solution of a symmetric generalization of equation ([2]) on the uniform lattice 
x{s) = s. Computational aspects of these new polynomials are described in detail 
giving their explicit representation as well as the three-term recurrence relation 
they satisfy. A full classification of weight functions and orthogonality supports is 
given together with computing the moments of the aforesaid weights. From this 
class all classical symmetric orthogonal polynomials of a discrete variable can be 
recovered (section |6|), and its limit relation with the continuous type of generalized 
classical symmetric orthogonal polynomials is given (see remark [2]). 

A regular Sturm-Liouville problem of continuous type is a boundary value prob- 
lem in the form 

(3) ^(k{x)^^^]+{Kg{x)-q{x))y„{x)=0 {k{x) > 0, g{x) > 0), 



dx \ dx ^ 

which is defined on an open interval (a, 6), and has the boundary conditions 

(4) aiy(a) + /3iy'(a) =0, a22/(fe) + /32y'(6) = 0, 

where ai,a2 and /3i,/32, are given constants and fc(a;), k'(x), q{x), and g{x) in ([3]) 
arc to be assumed continuous for x € [a,,b]. If one of the boundary points a and 
b is singular (i.e. fc(a) = or k(b) = 0), the problem is transformed to a singular 
Sturm-Liouville problem. 

Let Un and ym two eigenfunctions of the operator D[k[x)D) — q{x)I, where D 
is the standard derivative operator. According to Sturm-Liouville theory |29| . they 
are orthogonal with respect to the weight function q{x) under the given conditions 
(jll and satisfy the orthogonality relation 

g{x)yn{x)y„i{x)dx = / g{x)yl{x)dx Sn,„i- 



Many of special functions are orthogonal solutions of a regular or singular Sturm- 
Liouville problem having the symmetry property (<?!)„(— x) = (— l)"(/)„(a;)) so that 
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have found valuable applications in physics and engineering, as already mentioned. 
In [20], the classical equation ^ is symmetrically extended as follows. 

Theorem 1. [20j Let (j)n{—x) = (— l)"(/)„(x) be a sequence of symmetric functions 
satisfying the equation 

(5) A{x)(b':,{x) + B{xW„{x) + {\nC{x) + D{x) + anE{x)) 4>n{x) = 0, 
where 

1 — (—1)" Jo, n even, 



(6) 



1, n odd, 



and X„ is a sequence of constants. If A{x), {C{x) > 0), D{x) and E{x) are even 
functions and B{x) is odd then 



Q* {x)4)n{x)4)m{x)dx = [^J Q* {x)(t)„{x)dx j Sn,,n, 

where 

(7) = C(.)exp '^^^^ > dxj = -A-Zexp [^J ^dx 

The weight function defined in ^ must be positive and even on [—v, v\ and the 
function 

A{x)K{x) = A(x)exp i^j M^b-^da:) = exp (^j ^dx 

must vanish at x = v , i.e. A{v)K{v) = 0. In this way, since K{x) = g*{x)/C{x) 
is an even function so A{^i')K{—i') = automatically. 

Using this theorem, many symmetric special functions of continuous type have 
been generalized indnilSniinillllMlllIlllSl- 

Orthogonal functions of a discrete variable can similarly be solutions of a regular 
or singular Sturm-Liouville problem of discrete type in the form [15] 

(8) AiK*ix)Vy„{x)) + {X„wix)-q*ix))yr,ix)=0 > 0, u;(a:) > 0), 
where 

A/(a;)=V/(x+l) = /(x + l)-/(.T), 

and ([5]) satisfies a set of discrete boundary conditions like This means that if 
yn{x) and ym{x) are two eigenfunctions of difference equation (|S]), they are orthog- 
onal with respect to the weight function g* (x) on a discrete set [55] . 

Recently in [23] we have presented the following theorem by which one can gen- 
eralize usual Sturm-Liouville problems with symmetric solutions in discrete spaces. 
As a very important consequence of this theorem, we can introduce a basic class of 
symmetric orthogonal polynomials of a discrete variable with four free parameters. 

Theorem 2. [23] Let 4>n(~x) = (— l)"0„(a;) be a sequence of symmetric functions 
that satisfying the difference equation 

(9) A{x)l\V^n{x) + {A{-x) - A{x)) A(bn{x) 

+ (A„C(.t) + D{x) + <J„E{x)) c^„{x) = 0, 
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where, as usual, AV = A — V. If A(x) is a free real function and {C{x) > 0), D{x) 
and E{x) are even functions then 



p-i 

Y,W*{x)4>n{x)^.n{x) 



■p-l 

Y,w*{x)4>l{x) 



where 

(10) W*{x) = C{x)W{x) 

and W{x) is the solution of the Pearson difference equation 
A [A{x)W{x)] = [A(-.t) - A{x)] W{x), 

which is equivalent to 

W{x + 1) ^ A{-x) 
^ ' W{x) A{x + iy 

Moreover, the weight function defined in lilU\) must be even over one of the four 
following symmetric counter sets 

i) 5i = {—a — n, — a— n + a— 1,— a,a, a+l,...,a + n}, ogR, 

ii) 5*2 = 5*1 U {0} (as any odd function is equal to zero at x = 0), 

iii) ^3 = {. . . , — a — 71, — a — n + 1, . . . , — a — 1, — a, a, a + 1, . . . , a + n, . . . } (an 
infinite set) 

iv) 54 - ^3 U {0}, 

and the function A(x)W{x) must also vanish at x = a and x = 1 — jS , where 

e {5i, 52, 53, 54}. 

2. A BASIC CLASS OF SYMMETRIC ORTHOGONAL POLYNOMIALS OF A DISCRETE 

VARIABLE USING THEOREM [2] 

As a special case of equation ([5]), the following differential equation is defined in 

m- 

(12) 

x2(p*x2+g*X(x)+x(r*x2 + s*X(a;)-(n(r* + (n-l)p*)x2 + CT„s*)$„(x) = 0. 

Here one of the basic solutions is the symmetric class of orthogonal polynomials 
denoted by 



r s 

p q 



[n/2] 
fe=0 



[n/2] 
k 



'[ri/2]-(fc+l) 



n 

1=0 



(2^ + (-1)"+! +2[n/2])p* +r* 
{2i + (-1)"+! + 2)q* + s* 



-2k 



By referring to theorem [1] we observe in (fT2|) that A{x) = x^{p*x'^ + 9*) is a 
polynomial of degree at most four, C{x) = is a symmetric quadratic polynomial, 
D{x) = 0, and E{x) = s* . Since discrete orthogonal polynomials have a direct 
relationship with continuous polynomials, motivated by (jl2p . we suppose in the 
main equation ([9]) that 
4 

(13) A{x) = ^ aix\ C{x) = C2X^ + cq, D{x) = 0, E{x) = eo, 

i=0 



A BASIC CLASS OF DISCRETE ORTHOGONAL POLYNOMIALS 



5 



though there may exist some other possible cases. 

By noting the assumptions (|13p , we are now interested in obtaining a symmetric 
orthogonal polynomial solution. Hence, let 

(14) 0„(a;) =x" + <5„a;"-2 + ... , 
satisfy a three term recurrence relation as 

(15) 0„+i(a;) = x(l)n{x) - 7„(?!)„_i(a;), (with 0o(a;) = 1, 0i(x) = x). 
From ([9]), ([T3|) and ([M]). equating the coefficient in x"''"^ gives 

n (2a3 — a4{n — 1)) 



(16) 



C2 



provided that C2 ^ and also jasj + |a4| ^ 0. 

By using the eigenvalue A„ in (jf 6p and equating the coefficient in .t" we obtain 



- {6c2 (4ain - 2a2(n - l)n + ((-!)" - 1)) 

+04(71 — \)n (12co — C2{n — 3)(n — 2)) 
+Aazn (C2(n - 2)(n - 1) - 6co)} / {24c2 (04(8 - 2n) + 203)} . 
Also from and (|15p we have 

+ (5„+ia;"-i + • • • = X (x" + (5„a;"-2 + ...)- ^„ (x"-! + 5„_ia;"-3 + . . . ) ^ 
which implies 

(17) 7„ = (5„ - (5„+i. 

Therefore, in order that (pnix) is a solution of the following extra conditions 
must be considered for the initial values of n: 



eo = 2ai 



2a3Co 

C2 



Co ((2a3 - 04) Co + (02 - 2ai) C2) 



c2 



(18) 











+ fli {x^ 


4) 







As a conclusion, we get 

A{x) = a^x^ + 032;'^ + 02 l^a;'^ 

B{x) = A(-x) - A(x) = -2x (asx^ + ai) 
C{x) - CO (1 - 4a;2) , 
L»(a;) = 0, 

S(a;) = i (4ai +03), 



C2 = -4co. 



03 _ 04 
8 16' 



, (04(71 — 1) — 2a3) 
'^0 ' 



and 



(19) <5„ = {12ai (271 + (-1)" - 1) + 3a3 ((-1)" - 1) 

+71 (-12a2(?i - 1) + 203(2(71 - 3)71 + 7) 

-04(71 - l)((7i - 5)71 + 9))} / {2404(3 - 27i) + 4803} . 

For simplicity, if we set 

04 = 2a, 03 = o + 26, 02 = 6 + 2c, oi = c + 2d, 
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then pS]) changes to 



A{x) 
B{x) 

C{x) 



{2x + l){ax^ + bx^ + CX + d), 
-2x {{a + 2b)x'^ + c + 2d) , 



1 

4 

D{x) = 0, 

E{x) = i {a + 2b + 4c + 8d), 
Xn = 2n{an - 2{a + b)) for \a\ + \b\ ^ 0. 

Hence the following difference equation appears 

(20) (2x + \)(ax^ + bx^ + cx + d)AV0„(x) - 2x (a;^(a + 26) + c + 2d) A0„(x) 

1 o \ /a 



+ 2ri,(ari - 2(a + b)) 



a„ \- + b + 2c + 4:d 



0„(a;) = 0. 



Since the polynomial solution of equation (|20p is symmetric, we use the notation 
0„(x) = ^ " ^ 

for mathematical formulae and S'„(a, &, c, d; x) in the sequel. This means that from 
now we deal with just one characteristic vector V = (a, b, c, d) for any given sub- 
case. 

If the polynomial sequence {S'„(a, &, c, d; a;)}„>o satisfies a three term recurrence 
relation of type (fTSI) then by referring to (IT71) and we obtain 



Proposition 1. The coefficient 7„ m the three term recurrence relation U5\} is a 
rational expression in n where the numerator is a polynomial of degree 4 ciiT-d the 
denominator is a polynomial of degree 2, given by 



(21) 
where 



a b 



Y,i=oKi{a,b,c, d)j 



32(&-a(n-2))(6-a(n-l))' 



-Za , 



72 



and 



Ki{a, b, c, d) 
Kz{a,b,c,d) = 4a(3a + 25), 
K2{a, b, c, d) = -8 (3a2 + a(46 + c) + 6^) , 

i4:i(a, 6, c, d) = 2(3a + 26)(3a + 4(6 + c)) - 2a(-l)"(a + 26 + 4c + 8d), 
Koia, 6, c, d) = ((-1)" - 1) (3a + 26)(a + 26 + 4c + 8d). 

For n — 2m and n — 2m + 1, 7n in (|21ll becomes 

m (-a2(m - 1)'"^ + a (26(m - 1)^ + c(l - m) - d) + 6(6(1 - m) + c)) 



72m+l 



(6 - 2a{m - 1))(6 - a(2m - 1)) 

(a(m — 1) — 6) (— am'^ + 6m^ — cm + d) 
(2a m - (a + 6)) (2a TO - 6) 
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Remark 1. Once we have explicitly determined in the recurrence relation il5]}. 
a discussion about the situation of this coefficient is extremely important. For in- 
stance, analyzing the location of the zeros of orthogonal polynomials would give rise 
to a positive definite case when 7„ > (Vn G N), the quasi- definite case when 
7„ 7^ 0, and weak orthogonality case when 7„ = for some values of n. However, 
this discussion completely depends on the four parameters a, b, c and d, because 7„ 
is in general a rational expression in n. In section^ we provide factorized repre- 
sentations of "fn in two hypergeometric sequences, which would allow us to easily 
determine the orthogonality situation based on the parameters. 

Theorem 3. The explicit form of the polynomial S'„(a, 6, c, d; x) is 

ln/2] 



(22) 5„ 



a b 
c d 

'["/2]-l 

n 



x-^(^l) 

J=0 



[n/2] 

i 



a{i + - Hi + '^n)^ + c{i + ct„) - d 
a{i + [n/2] - ct„+i) - b 



(cFn - x)j{an + x) j , 



where [x] denotes the integer part of x, (A)„ ~ A{A + 1) • • • (A + n — 1) = T(A + 

-1 

n)/r{A) forn>l with {A)o = I and JJ (.) = 1. 

4=0 

Moreover, since {—x)j{x)j = (—1)-' rife=o(^^ ^ ^'^)' ^ ^""^ ^ -^"^ + ^ 

(fH) becomes 



2m 



a b 
c d 



and 

S2m-\-l 



a b 
c d 



= p^ix')^Yl 

3=0 



X ) = xRm{x^) 



n 



ai^ — bi^ -\- ci — d 
a{i + TO — 1) — 5 




? 7j n 



a{i + 1)3 - b{i + 1)2 + c{i + l)-d 



1=3 



aii + to) 



\{{x'-e). 



k=l 



Proof. By defining the symmetric basis 



(23) K{x) = (-l)!"/^]^.-. (^^ _ + x)[„/2] 

[n/2]-l 



fc=0 



the following straightforward properties are derived 

(24) AMx) = n^„_i(x) + !fcll^„_2(a;) + ~ ^ ^„_3(x), 

(25) x??„(x) i?„+i(a-) + a-„+i (^-j ?9„_i(2;), 

(26) AVi9„(x) = n(n - l)K-2{x). 
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The expansion of 5„(a, b, c, d; x) in terms of the above basis yields 

[„/2] 

(27) 5™ (a, 6, c, d; x) = ^ Cj {n)-d2j+cr„ (x) 
j=o 

' m 

^{-iycj{2m){-x)j{x)j, n = 2m, 

- < 

\ m 

n = 2m + 1, 

then by using the Navima algorithm [5J [301 [31] we can reach a solvable recurrence 
relation for the connection coefficients Cj(n). So, by using the latter properties 
(IMl)~(I2ni) and substituting (gT]) in ^ we obtain 

[n/2] 

= ^ Cj{n) {Ej{n) + i?2j+<t„+2 + Fj{n)d2j+a-„ + G.j{n)-d2j+a„-2) , 
j=o 

where 

Ej [n) = 2 {2j + cr„ - n) (a(2j + n - 2) + a(T„ - 25) , 

F,in)=4aij-lfj{2j-l)cTn-i 

+ i(T„+2 {2j + CT„ + 1) 2 (2j + (T„ - ?i) (a(2j + n - 2) + aCT„ - 26) 
+ 2ja„+i (a (4/ _ Gj^ - j{n - 3)(n + 1) - l) + 2fe(j(-4j + n + 3) - 1)) 
+ I i<^n (ffn (cr„ (a„ (8aj + acr„ - 4a - 2fe) + a(24(j - l)j + 5) + 46(1 - 3j)) 
+2a(2j - l)(8(j - 1).? + 1) + 86(2 - + Ac) + a (4(j - l)j(l - 2jf + l) 
+4j(6(l - 4(j - + 4c) - 4c) - 2n(a + 6) + an^ 

-4j(6 + 4(c + d)) + 8,f{b + 2c)) , 

and 

Gj{n) = U - l)i(2j - 1) (2(i - l)a„_i (2a(j - l)V„_i + 6 + 2c) 
-2a„+i ((j - l)2(a + 26)(T„_i + c + 2rf)) 

+ I (2j + a„ - 1) i2j + a„) (a„ (2j + a„ - 1) (a„ (2 ((j - 1) (-4j(a + 6) 

+4aj^ + a) + 2c) + dn ((t„ {6aj + a(T„ — 4a — 26) 
+ {2j - l)(a(6j - 5) - 46))) + 6(4j - 2) - 8(c + d) + 8cj) + 8d) . 

But since z?n(a;) is linearly independent, the coefficients Cj{n) would satisfy the 
relation 

Ej-i{n)cj^i{n) +Fj{n)cj{n) + Gj+i{n)cj+i{n) = 0, 

which is explicitly solvable with the initial conditions Cmin) = for m > [n/2], and 
C[n/2](?^) = 1, providing ([22]). □ 
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Remark 2. By substituting the characteristic vector 



(a, b, c,d) = ( ^h, — — 



' +r* 1 g* + 5" 

and X ~ t/h in relation i22\} . the corresponding polynomial satisfies a difference 
equation of type I12U\} that formally tends to the differential equation U2\) as /i — > 0. 
Hence, the following limit relation can be directly deduced: 



lim Sn 



/ i 
V 



p 



q* 2q* 
1 q* + s* 

h 



2q*h 




As the recurrence relation (jlSp is explicitly known, the complete form of the 
orthogonality relation is 



(28) J2 ^ 



a b 
c d 



a b 
c d 




a b 
c d 

a 
c 



Sn- 



where 



W 



a b 
c d 



W{x), 



is the original weight function and W(x) satisfies the difference equation 



(29) 



W{x + 1) (1/2) -X 



-ax 



bx^ 



cx 



W{x) (3/2) + xa{x + 1)3 + b{x + 1)2 + c{c+l) + d' 



By noting that A{x) = {2x + l)(o 



d) for \a\ + |6| 7^ 0, two cases 



can generally happen for the parameter a in (j29p . i.e. when a ^ and b arbitrary 
or a = and 6 7^ 0. 

In the first case, since any arbitrary polynomial of degree 3 has at least one real 
root, say x = p e R, the aforementioned A(x) can be decomposed in three different 
forms, i.e. 

(30) A{x) = {2x + 1) {x -p){ax'^ +UX + i) 

{2x + 1) (x — p) (a x2 + u x + v) , (w2 < 4au), 
{2x + 1) a{x ~ p){x ~ q){x — r), {u^ > 4av), 
(2a; + 1) a(x -p)(x - q)2, (^2 = 4ati). 

Similarly, in the second case when a = and 7^ 0, A{x) can be decomposed as 

(31) A{x) = {2x + l) {bx^ +cx + d) 

' {2x + l) {bx'^ +cx + d) , (c2<4M), 



{2x + l)b{x-p){x-q), 
^{2x + l)b{x-p)^, 



(c2 > 4bd), 
(c2 = Abd). 



For the two sub-cases A{x) = (2.T + 1) (x—p) (ax 



v) in dSni) and A{x) = 



{2x + 1) (6a;2 +cx + d) in (j3T|) . the difference equations corresponding to ([29|) 
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respectively take the forms 

(■xo\ W{x + l) ^ {l/2)-x p + x ax^-ux + v (2 . ^ 

W{x) {3/2)+x p-x-la{x + l)^ + u{x + l) + v 

and 

W{x + l) _ {ll2)-x bx^-cx + d 

^ ' W{x) {3/2) + x bix + l)^ + c{x + l) + d 

Since the denominators of the two fractions ([5^ and l\'S'3\i are not decomposable 
in R, we shall deal with them in a separate work [under preparation]. However, 
it is important to note that the cases analyzed in this paper allow us to recover 
all classical symmetric orthogonal polynomials of a discrete variable (section ^ . 
Hence, let us consider the rest of cases. For the second sub-case of without 
loss of generality take a = 1 and then A{x) = {2x + l)(x — p){x — q){x — r) where 
p^q^r € R, which indeed covers the third sub-case too. For the second sub-case 
of (|3ip we can similarly consider Aix) = {2x -I- l)(a; — p){x — q) where p,q £ R. 
This means that there exist only two orthogonal sequences of Snia, b, c, d; x) when 
A{x) is decomposable. Moreover, there exist 20 different weight functions on 32 
supports for these two sequences, as the weight functions are not in general unique 
leading to an indeterminate moment problem. In summary, when the polynomial 
A[x) is of degree four, a three parametric family (with 14 subcases) appears and 
when A{x) is of degree three, a two parametric family (with 6 subcases) would 
appear. Finally, when A{x) is of degree 2, the well known classical symmetric 
discrete families appear. 



3. Orthogonality supports 

The question is now how to determine the restrictions of the parameter 9 in the 
orthogonality support [—6', 0] in ((28|) ? To answer, we should reconsider the main 
difference equation on [a, /? — 1] = [—6', 6] and write it in a self-adjoint form to 
eventually obtain 

(34) A (^(x)M^(x) ((/)„, (x)V0„(x)-</.„(x)V0„ (x))) 

x=-e 

+ (A„-A„) ^ C^-xAw{x)(t>n{x)(t)m{x) 

x=-e ^ ^ 
(-1)" - (-1)" (a 



+ 5 + 2c + 4d) W{x)c^n{x)4>^{x)^Q. 



2 \2 

x = -U 

On the other side, the identity 

(j)„i{x)\'(f>nix) - (i)n{x)V 4>,n{x) = (/)„ (a;)(/),„ (x - 1) - (l>m{x)(j)„{x - 1), 
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simplifies the first sum of ([34|) as 



e 



(35) A(^(a:)M^(a:)(0™(x)V(/)„(x)-0„(a:)V(/),„(a:))) 




= A{x)W{x)[(j)n{x)(t)rn{x-l)- (/),„ (a;)(^„ (x - 1) ) 



= A{e + i)w{e + 1) {<Pn{e + 1)</)™(0) - + 1)<^„(0)) 

- A{-e)w{-e) iM-9)<Prni-e - 1) - 0„^(-^^)</)„(-0 - 1)) . 



By taking into account tfiat all weight functions are even, i.e. W{—x) = W{x), 
the polynomials are symmetric, i.e. 4>nix) = (— cc), and the Pearson differ- 
ence equation (HH) is also valid for x = 9, i.e. A{e + 1)W{0 + 1) = A{-9)W{0), 
relation ([35]) would be finally simplified as 



Since (j)ra{0)4>n{d + 1) ^ 4'n{9)4>m{9 + 1) 7^ 0, two cases can in general happen for 
the right hand side of (p6)) : 

i) If ri + m is odd then 1 + (— 1)'^+'" = q ^j^j (ggj) is automatically zero. 
However, this case is clear as the sum of any odd summand on a symmetric 
counter set is equal to zero. 

ii) If simultaneously A{—9) = and W{9) ^ 0, then ([36]) is again equal to 
zero. This condition is in fact the key point for finding the orthogonality 
supports. 

4. Two HYPERGEOMETRIC ORTHOGONAL SEQUENCES OF S'„ (a, 6, C, d; x) 

In this section, we introduce two hypergeometric sequences of symmetric orthog- 
onal polynomials, which are particular cases of S'n(a, 5, c, d\ x) corresponding to two 
aforementioned cases a 7^ and a = 0, respectively and then obtain all possible 
weight functions together with orthogonality supports for these two sequences. 

4.1. First sequence. If the characteristic vector 



6 



(36) A(^(a;)VF(x)(0„,(x)V</)„(x)-(/.„(x)V(/)„(.T))) 



= A{-e)W{9) (1 + (-1)"+'") (0„.(0)0„(0 + 1) - <^„(0)<^™(0 + 1)) . 



(a,fe, c. 



d) = {l,-{v + q + r),pq + pr + qr,-pqr), 



for p,q,r E R. is replaced in (PO]) . then 



{2x + l)(a; - p){x - q){x - r)AV0„(2:) 

— 2x (a;^(l — 2p — 2(7 — 2r) + pq + pr + qr — 2pqr) A(/)„(a;) 
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has a polynomial solution, which can be represented in terms of hypergeometric 
series as 



(37) Sn 



1 —{p + q + r) 

pq + pr + qr —pqr 

_{p + Crri)[n/2](<? + CT„)[„/2](7' + Cr„)[„/2] 



X x^-aFz 



([n/2] + P + g + r- l + (7„)f„/2] 
-[n/2], [n/2] +p + 9 + r - 1 + o-„, cr„ - x, ct„ + a; 



where (t„ is defined in (|6]) and 4-^3 is the well known hypergeometric function of 
order (4,3). Moreover, it satisfies the recurrence relation 



(38) (j>n+l{x) = X(j)n{x) ~ 7„ 



1 —{p + q + r) 

pq + pr + qr —pqr 

{(j>o{x) = l, (t)i{x) = x), 



where 



(39) 7r 



{-2n* - 4(-3 + 2p + 2q + 2r)r 



1 —{p + q + r) 

\pq + pr + qr —pqr 

-8 {p^ + p{3q + 3r - 4) + 3qr + {q - 4)g + - 4r + 3) 

+ (2(-l)"(-l + 2p){-l + 2q){-l + 2r) 

-2{-3 + 2p + 2q + 2r)(3 + 4g(-l + r) - 4r + 4p(-l + q + r jj) n 

+ {-1 + (-!)")(-! + 2p){-l + 2q){-l + 2r)(-3 + 2p + 2q + 2r)} 

I {32(n + p + q + r -2){n+p + q + r -\)} , 

dividing to 

n{n + p + q — \){n + p + r — \){n + q + r — 1) 



and 



72n = 



72n+l 



{2n + p + q + r -2){2n + p + q + r -\) 
{n + p){n + q){n + r){n+p + q + r -\) 



{2n + p + q + r- 1){272 + p + q + r) 
The latter representations would allow us to analyze the sign of 7,1 in terms of the 
values of p, q and r as explained in Remark [T] 

By noting the relations (|38p and (|39p . the orthogonality relation of the first 
sequence is 



1 —{p + q + r) 

pq + pr + qr —pqr 



^ I Sn{x)Sm{x) 



n 

fc=i 



7fe 



1 —{p + q + r) 

pq + pr + qr —pqr 



1 -{p + 1 + r) 

pq + pr + qr —pqr 



in which 



Sn {x) — Sn 



1 —{p + q + r) 

pq + pr + qr —pqr 
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and, finally 



W 



1 -{p + q + r) 

pq + pr + qr —pqr 



denotes the original weight function corresponding to the hypergeometric poly- 
noniialdenotes the original weight function corresponding to the hypergeometric 
polynomial (|37p . As a consequence, we obtain 

Proposition 2. The function W{x) satisfies a particular case of the difference 
equation Ii29\) 

, . W^(a; + 1) _ -a: + (1/2) -x - p -x - q -x - r 

^ ^ W{x) ^ x + (3/2) x + l-px + 1- qx + l-r' 

giving rise to the 16 symmetric solutions of listed below 

(41) Wi{x-p,q,r) = {T{l-p + x)T{l-p-x)T{l-q + x)T{l-q~x) 

xr(l - r + .T)r(l - r - a;)r(3/2 + x)T{i/2 - x))~^ , 



(42) 


W2,i[x;p,q,r) 












x)V{p — x) 








T(l-q + x)T(l -q- x)T{l - r 


+ x)T{l-r -x)T{i/2A 


-a;)r(3/2- 


xY 


(43) 


W2,2 {x; p, q, r) = W2,i{x; q, p, r) 










r(g + 


x)V{q - x) 








r(l -p + x)T{l -p-x)T{l~r 


+ x)r(l -r-a;)r(3/2H 


-a;)r(3/2- 


xY 


(44) 


W2,3 {x;p, q, r) = 14^2,1 (2:; r, q, p) 










r(r + 


x)V{r — x) 







r(l - p + a;)r(l -p- x)T{l -q + x)T{l -q- a:)r(3/2 + a;)r(3/2 - x) '' 
T{p + x)T{p — x)T{q + x)r(q — x) 



(45) W3,i(s;p,9,r) 
(46) 



W3,2ix; p,q,r) = W3^i{x;p,r,q) 
(47) 



r(l - r + a;)r(l -r- x)T{3/2 + x)T{3/2 - x) ' 
T{p + x)T(j> — x)T{r + x)T(r — x) 



W3,3{x]p,q,r) = W3^2{x\r,q,p) 
(48) Wi,i[x-p,q,r) = 



T{l-q + x)T{l ~q- x)T{i/2 + x)r(3/2 - x) ' 

V{q + a;)r(g - x)T{r + x)T{r - x) 
T{l-p + x)V{l -p- x)V{Z/2 + a;)r(3/2 - x) ' 

r(p + x)T{p - a;)r(-l/2 + a-)r(-l/2 - x) 



r(l - «7 + x)T{l -q- x)T{l - r + x)T{l - r - x)'' 



(49) Wi.2{x]p,q,r) = Wi.i{x;q,p,r) 

T{q + x)T{q-x)T{-l/2 + x)T{-l/2^x) 



T{l-p + x)T{l -p- x)r{l -r + x)r{l -r-x)' 
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(50) W4^3ix;p,q,r) = Wi,i{x;r,q,p) 

r(r + x)r{r - x)T{-l/2 + x)r{-l/2 - x) 



(51) 

W5,i{x;p,q,r) = 



T{l-p + x 


)T{l-p-x)Til 


- q - 


^x)T{l 


- q 


-x) 


r{p + x)T{p — x)r{q 


+ x)T{q-x)T{~ 


1/2 4 


-x)r(- 


1/2 


-x) 


r(i 


— r + 2^)r(l — r 


-x) 








) = W5^i{x;p,r,q) 
T{p + x)T{p — x)T{r 


+ x)T{r - x)r{- 


1/2 H 


^x)T{- 


1/2 


-x) 


r(i 


-q + x)ril-q 


-x) 








) = W5,i{x;r,q,p) 
r(g + x)r{q - x)T{r 


+ x)r{r - x)r{- 


1/2 ^ 


-x)T{- 


1/2 


-x) 



(54) Wf,{x-p,q,r) 



T{1 - p + x)T{l - p - x) 
r(-l/2 + x)r(-l/2-x) 



r(l - p + a;)r(l - p- x)T{l -q + x)T{l -q- x)T{l - r + x)T{l - r - x)'' 
T{p + x)T(j) — x)T{q + x)T{q — x)T{r + x)T{r — x) 



(55) W7{x;p,q,r) 



r(3/2 + j:)r(3/2-.T) 



(56) W8{x;p, q, r) = T{p + x)T{p - x)r{q + x)r{q — x) 

X r(r + x)r{r - a;)r(-l/2 + x)r(-l/2 - x). 

As the original weight functions corresponding to aU above 16 cases are as 

^-x'^Wkix), 
the two following identities arc remarkable in this direction 

(1/4 - x^)r{-i/2 + x-)r(-i/2 -x)^ r(i/2 + x)r(i/2 - x), 

and 

1/4 _ 1 
r(3/2 + a;)r(3/2 - x) ~ r(l/2 + x)r(l/2 - x) ' 
For example, for the last given case the original weight function becomes 



1 -{p + q + r) 

pq + pr + qr —pqr 



X \ = \ x^ ] W^{x;p,q,r) 



= T{p + x)T{p - x)T{q + x)T{q - x)T{r + x)T{r - x)V \^+xjT\^-x 

By noting section 3, the following table now shows the orthogonality supports of 
each weight function given in (|4ip to (|56p together with their parameter restrictions 
in which Z = {0,-1,-2,...}. 
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Table 1 . Orthogonality supports for the first hypcrgcomctric sequence 



Wk{x) Support Parameter restrictions 



Wi{x;p,q,r) 




pe 


z 


, 1 - g±p ^ 


Z 


, 1 - r ±p <^ 


Z . 








z 




z 


1- r ±q <^ 


Z . 




[— r, r] 


r G 




, 1 - q ± r ^ 


Zj 


, I - p±r ^ 


Zj 


W2,i{x;p,q,r) 






z 


,p±q^Z- 


, 1 


~r±q(^Z- 






[— r, r] 


r e 


z 


,p±r<^Z- 


, 1 


-q±r ^Z- 




W3^i{x;p,q,r) 


[— r, r] 


r e 


z 


,p±r^Z- 


, q 


±r (^Z~. 




W4,i{x;p,q,r) 




q& 


z- 


,p±q^Z- 


, 1 


-r±q^Z- 






[-r, r] 


r e 


z- 


,p±r^Z- 


, 1 


~q±r ^Z- 




W5^i{x;p,q,r) 


hr, r] 


r e 


z 


,p±r(^Z- 


, q 


±r (^Z-. 




W6{x;p,q,r) 




pe 


z 


,l-q±p<^ 


z 


-,l-r±q<^ 


z . 




[-9. <?] 


q € 


z 


,l-r±q(^ 


z 


-,l-p±q<^ 


z . 




[-r, r] 


r e 


z- 


, 1 — q±r ^ 


z 


~ , 1 — p±r ^ 


z . 


W7{x;p,q, r) 
















W8{x;p,q,r) 

















Remark 3. Since some weight functions are symmetric with respect to the pa- 
rameters p, q and r, e.g. W4,^2ix;p,q,r) = W4,^i{x; q,p,r) and W4,^3{x;p,q,r) = 
W4^i{x]r,q,p), their orthogonality supports and parameter restrictions can be di- 
rectly derived via the rows of table{^ by just interchanging the parameters. Also note 
that Wj{x;p,q.,r) and Ws{x;p,q,r) have no valid orthogonality support. Therefore, 
there are totally 24 eligible orthogonality supports for the first sequence. 

Remark 4. Let the weight functions corresponding to the first sequence be indicated 
as Qi{x) ~ (1/4— x'^)Wi{x). Then they satisfy the difference equation 

Qi{x + l) _ {p + x){q + x){r -\- x) 

Qi{x) ^ {-p + x + l){-q + x + l){-r + x + iy 

which is equivalent to 

A{Mi{x)g,{x)) = Ni{x)g,{x), 

where 

Mi{x) ^ {x - p){x - q){x - r), Ni{x) = 2x'^{p + q + r) + 2pqr. 

Thus, the hyper geometric polynomials \31\ l constitute a A-semiclassical sequence of 
class one [18] . for which we have explicitly given their orthogonality condition, three 
term recurrence relation, hypergeometric representation and in the sequel in section 
O we will compute the factorial moments with respect to the basis iSn {x) defined in 

mi. 
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4.1.1. A numerical example for the first sequence. Let us replace p = —9, q = —10 
and r = — 11 in (l37l) to set 



(57) Sn 



( 1 


30 




299 


990 






(-9 + cr„)[„/2](-10 + cr„)[„/2](-ll + cr„)[„/2] 



X X"" 4^3 



(-31+ [n/2]+a„)[„/2] 
-[n/2], -31 + [n/2] + o-„, cr„ - x, a„ + x 
-9 + cr„, -10 + (T„, -ll + cr„ 



Since Z^, 1 — qip^ and 1 — r±p^Z^, referring to Table [T] shows 
that there are two possible weight functions for the selected parameters which 
are orthogonal with respect to the sequence S'„(l, 30, 299, 990; x) on the support 
{-9, -8,..., 8, 9}, i.e. 



Wi 



( 1 


30 






299 


990 


X 





x^j l¥i(a;;-9,-10,-ll) 

= (r(io + x)r(10 - a;)r(ll + a;)r(ll - x) 

r(12 + a;)r(12 - a;)r(l/2 + a;)r(l/2 - a;))"^ 



and 



We 



( 1 


30 






299 


990 


X 





r(l/2 + a;)r(l/2-a:) 



r(io + a;)r(io - j;)r(ii + a:)r(ii - x)r(i2 + x)r(i2 - x) ■ 

Hence there are two orthogonality relations corresponding to the polynomials (|57p , 
respectively as follows 



x=-9 



( 1 


30 




299 


990 






X S'„(l,30,299,990;a;)5„(l,30,299,990;x) 



36183421612800000(-l)"(-29) _^ (-19) ^ „ j _^ (-18) ^ „ j _^ (-17) ^ „ j _^ 



(-31)„(-30)„ 
x(-9)L^j_,(-8)L^j_,(-7)L^j_,r( 

where [x\ denotes the floor function, and 



+ l]Sn 



a:=-9 



( 1 


30 






299 


990 


X 





13/tt, 1 = 1, 
f3n, 1 = 6, 



in which 



/3 



667 



1998530094928466929986605067918114816000000000 
4.2. Second sequence. If the characteristic vector 

{a,b,c,d) = {Q,l,-p- q,pq), 
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for p, q e R, is replaced in the difference equation (|20l). then 

{2x + l)ix - p){x - q)A\'(f>nix) - 2x (2x^ + p{2q - 1) - q) A(f>„ {x) 



Ani - - x' 



(2p-l)(2g-l)U„(x)=0 



.4 7' 2 

has a basis solution, which can be written in terms of hypergeometric series 



(58) 5„ 



1 

-p-q pq 



X ^ = (p + cr„)[„/2]((7 + cr„)[„/2] 



X x^" 3F2 



Ir 



q = \ (-2n' - 4n(p + q-l) + {{-IT - 1) [2p - f )(2g - 1)) 



[n/2], cr„ - X, a„ + x 

P+(Tn, q + (Jn 

The above polynomial is a limit case of the polynomial p7p when r — > 00. Moreover, 
it satisfies a recurrence relation of type (|15p with 

1\ 1 
which implies 

72n = + n+p + q), 

and 

72«+i = -(n + p){n + q). 
Hence, the orthogonality relation corresponding to the second sequence takes the 
form 



1 



^-^ \ \ —p — q pq 



X Sn 



1 

-p-q pq 



X 



S71 



1 

-p-q pq 



1 

-p-q pq 



1 



W 



\ — p — g pq 



Wix), 



X Sn- 











X 









fc=i 

where 

1 

-p-q pq 

is the original weight function. Thus, we obtain 

Proposition 3. The function W{x) satisfies a particular case of the difference 
equation i29\) as 

, , W{x + 1) _ 1/2 - X -x-p -X - q 

^ ' W{x) ~ 3/2 + xx + l-px + l-q' 

giving rise to 8 symmetric solutions for equation \59fl respectively as follows: 

(60) Wc>{x-p,q) 

1 

^ r(l - p + a;)r(l -p- x)V{l -q + x)T{l -q- x)T{S/2 + x)T{3/2 - x) ' 

r{p + x)T{p — x) 



(61) WioA^;P,<l)^ 



T{l-q + x)T{l -q- x)T{3/2 + x)T{3/2 - x) ' 
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Tiq + x)T{q~x) 



(62) WioAx;p,q) 



(63) Wi,ix;p,q) 

(64) Wi2{x;p,q) = 



r(l - p + a;)r(l - p - x)T{'S/2 + x)T{3/2 - x) ' 
r{p + x)T{p — x)r{q + x)T{q — a;) 



r(3/2 + x)r(3/2-x) 
r(-l/2 + a:)r(-l/2-a;) 



(65) W,sAx;p,q) 

(66) W,3Ax;p,q) 



r{l-p + x)r{l -p- x)r{l -q + x)T{1 -q-x)' 
V{p + x)T{p - x)T{^l/2 + x)T{-l/2 - x) 



V{l-q + x)V{l-q-x) 
T{q + x)T{q - x)T{~\/2 + x)r(-l/2 - x) 



T{l^p + x)T{l-p-x) 

(67) Wu{x-p, q) = T{p + x)T{p - x)T{q + x)T{q - x)T{-l/2 + x)T{-l/2 - x). 

Tabic [5] shows the orthogonahty supports of each above-mentioned weights and 
their parameter restrictions. 

Table 2. Orthogonahty supports for the second hypergeometric sequence 



Wk{x) Support Parameter restrictions 



Ws){x;p,q) 


[-P,P] 


pez- 


, 1 


-q±p<^ 


'.Z-. 






qez- 


, 1 


-p±q^ 


[ z . 


WwAx\p,q) 




qeZ- 


, p 


±q^Z- 




Wii{x;p, q) 












Wi2{x;p, q) 




p€Z- 


, 1 


-q±p<^ 


'.Z-. 




hi, q] 


qeZ- 


, 1 


-p±qf 


[ z . 


Wi3,i{x;p,q) 




qeZ- 


, p 


±q(^Z- 




Wi4{x;p, q) 













Remark 5. As the main orthogonality relation h28\) shows, an important part that 
one has to compute in norm square value is W(x; a, 6, c, d) where [—9,9] 

is the same orthogonality supports as determined in tables [I] and Since the 
functions {Wk{x; a, b, c, d)}1^^ are all even, the aforesaid sums can be simplified on 
its orthogonality supports by using two identities 

T{p + x)^Tip)ip), and Tip-x) = ' , 

i-P)x 
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and this fact that 
9 



^ W{x; a, &, c, d) = 2 ^ W{x] a, b, c, d) - W{Q] a, 6, c, d). 



For example, for the first given weight function Wi{x]p,q^r) on e.g. [—p,p\ we 
have 



(^i Wi{x;p,q,r) 



^r2(l - p)r2(l - q)T^{\ - r) \^ (1 - pUI - qUl - r). 



- 1 



Remark 6. Let the weight functions corresponding to the second sequence be indi- 
cated as Qj{x) = (1/4 — x'^)Wj{x). Then they satisfy the difference equation 



Qjix) 



{p + x){q + x) 



{~p + x + l){-q + x + l) 



which is equivalent to 



A{M2g,ix)) = N2gjix), 



where 



A'hix) ^{x- p){x - q), N2{x) = -2 {pq + x^) . 



4.2.1. A numerical example for the second sequence. Let us replace p = 14 and 
q = -10 in (1371) to get 



1 

-4 -140 



X ^ = (14 + cr,i)[„/2](-10 + cr„)[n/2] 



[n/2], On - X, an+ X 
U + a^, -10 + cr„ 



1 



Since q E Z~ and p ± q ^ Z~, referring to Table [2] shows that there are two 
possible weight functions for the selected parameters which are orthogonal with 
respect to the sequence §„(0, 1, — 4, — 140; x) on the support { — 10, — 9, . . . , 9, 10}, 
i.e. 



1^10,1 







1 

-140 



1 



x' W^io,i(a;;14,-10) 



(13 - a;)(12 - a;)(ll - x){x + ll){x + 12){x + 13) cos(7ra;) 



and 







1 



= (i-x2)W^i3,i(x;14,-10) 



-4 -140 

= 7r(13 - a;)(12 - x)(ll - x){x + ll){x + 12){x + 13) sec(7rx). 
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Hence there are two orthogonality relations corresponding to the polynomials 
([68| respectively as follows 



10 



a;=-10 



1 

-4 -140 



X j ^„(0, 1, -4, -140; x)5„(0, 1, -4,-140; x) 

i)n+i(„9)^^^r ([H^J + 15) r ([tJ + 1) r ([|J + 4) 

3736212480 



where 



'01 



a:=-10 



-4 -140 



fl0296/7r, j = 10, 
[I02967r, j^li. 



5. Moments of the two introduced sequences 

To compute the moments of a continuous distribution, different bases should be 
considered. For example, in the normal distribution the canonical basis {x-'lj^o is 
used to get 

g-xV2 fO> ri = 2TO + l, 

x'' — T^dx = < 2"+^/2 / 1 



'2-K , r m H — , n ~ 2m, 

while for the Jacobi weight function (1 — x)" {1 + x)^ as the shifted beta distribution 
on [—1, 1], using one of the two bases {(1 — xy}j or {(1 + xy}j is appropriate for 
this purpose. 

This matter similarly holds for the moments of discrete orthogonal polynomials. 
For instance, in the negative hypergcometric distribution corresponding to Hahn 
polynomials, it is more convenient to use the Pochhammer basis {(— a;)„}„>o, in- 
stead of the canonical basis, to get 

^ r{N)r{a + p + 2)r(a + N -x)T{p + x + l) 

^ r(a + i)r(/3 + i)r(a + i3 + n + i)r{N - x)r{x + 1) ^ 

(l-jV)„(/3 + l)„ 
^ ' (a + /3 + 2)„ 

Following this approach, for the weight functions gi{x) appeared in the two 
introduced liypergeometric sequences, we shall compute the moments of the form 

= E ^n{x)g^{x) = ^ dn{x) Q ~ xA W^{x). 

x=-e x=-e ^ ^ 

where the basis dn{x) is defined in (p3)) . 

Since '&2n+i{x) are odd polynomials, clearly all odd moments with respect to 
this basis are zero. Moreover, from (P7)) and ([55]) and also using the orthogonality 
property of the polynomials it can be proved by induction that the even moments 
corresponding to the first and second sequences respectively satisfy the following 
recurrence relations 

{p + n-l){q + n-l){r + n-l) 

[Qi)2n — ; ; ; z. l£'ij2n-2, 

p+q+r+n— 1 
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x=-0 



and 

{Q])2n = -(n + p - l)(n + q- l){Qj)2n~2- 

Therefore, if the aforesaid weight functions are normahzed with the first moment 
equal to one, then we eventuaUy obtain 

rO, n = 2m + l, 

x=-e ^ ^ — ■ ■ — r , n = 2m, 

K. {p + q + r)m 

for the weight functions of the first sequence, and 

\ fO, n = 2m + l, 

^ ' .(-l)™b)™(g)m, n = 2m. 

for the weight functions of the second sequence. 

6. A PARTICULAR EXAMPLE OF Sn{a,b, C, d; x) GENERATING ALL CLASSICAL 
SYMMETRIC ORTHOGONAL POLYNOMIALS OF A DISCRETE VARIABLE 

In this section, we present a particular interesting example of S'„(a, 6, c, d; x) 
that generates all classical symmetric orthogonal polynomials of a discrete variable 
studied and analyzed in [5] and is different from the two introduced hypergeometric 
sequences. 

If a = —2(6 + 2c + 4d) in the main difference equation ([20]). after simplification 
of a common factor we get 

{x^{b + 2c + Ad) + x{c + 2d) + d) AV0„(x) - 22;(c + 2d)A(t>„ {x) 

+ n{-{b{n - 1) + 2{n - 2)(c + 2d)))<j)„{x) = 0, 

which is the same equation as analyzed in [3], 

(ax^ + bx + c)AV2/„(x) - 2bxAy{x) + n(a(l - n) + 2b)y{x) = 0, 

for 

d = c, c = 6 — 2c, b ~ a — 2b, a = —2a. 
Hence, the particular polynomial 

[n/2] 

(69) 5„ 



-2a a -26 
6- 2c c 



X " 



"/2] /r ,„1 



'[n/2]-i (2^ + 2a„ + 1) (-a(^ + a„)2 + + _ a) \ 

n A \ ("^^ ^ 2:)j(cr„ + a;)j, 

i=i 2 f 6 + a(o-„+i - i - [7772])) - a J 

generates all cases of classical symmetric orthogonal polynomials of a discrete vari- 
able respectively as follows: 

Case 1. If a = 0, 6 = 1 and c free in (j69| then a{x) = {dx^ + bx + c) — x + c has 
one real root and the symmetric Kravchuk polynomials [5l [12] are therefore 
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derived as 



S„ 



-2 
1 - 2c c 



-n, — c — a:; 
-2c 



a; j =2-"(-2c)„2Fi ^ 

= fcii/2)(a; + c; 2c), 

which are orthogonal with respect to the weight function 
1 



g{x) 



for X g {— c, — c + 1, . . . , c — 1, c}, 



r(c-.T + l)r(c + a; + 1) 
when 2c G N. 

Case 2. If a = 1, 6 free and c = 6^/4 in ([M]), then cr{x) = (aa;2 + Sa; + c) = {x + b/2)'^ 
has a double real root and the symmetric Hahn-Ebcrlein polynomials [28j 



-2 1-26 
6(2-6)/2 P/4: 



(n- 26- 1) 
-n, 71 -26-1, -a; -6/2 
-6,-6 



1 



/i^o)(a- + 6/2,6+l). 



are orthogonal with respect to the weight function 

1 



q{x) 



r^{x + 1 + 6/2)r2(-2; + 1 + 6/2) 

for X € {-6/2, -6/2 + 1, . . . , 6/2 - 1, 6/2} when 6 e N. 
Case 3. If a = 1, 6 = —5i — 82 and c = 5x52 in (|69)) . then cr(x) has two different 
real roots 5i and 82 and depending on the values 5i and 82, the symmetric 
Hahn-Eberlein polynomials [31 [5S] are derived as 



Sn 



-2 l + 2(5i+52) 

-(52-(5i(l + 2,52) 5i52 



■ 3F2 



-n, n + 2(^1 + ^2) — 1, (^2 - a; 
5i + ^2, 2(52 



^ {6l+62)n{252)n 

{n + 2(,5i + ^2) - 1), 

^^(<52-5i,52-A-i)(2.„^2;l-252), 

when —2(52 <= N or the symmetric Hahn polynomials [3l[28] as 

-2 l + 2((5i+(52) 

-(52 - (5i(l + 2^2) 8182 

{5l+52)n (2(5i), 



(n + 2((5i+(52)-l)„ 



3-P2 



— n, n + 2(^1 + 82) — 1, (5i — a; 
<5i+<52,2(5i 
^ /i(fi+^2-i,5i+5i-i)(^ - <5i; 1 - 2,5i), 

when — 2(5i £ N. Finally, for (5i + (52 = 1 we get to Gram polynomials |13) . 
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